Dimensional Crossover in Quantum Antiferromagnets 
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The dimensional crossover in a spin-S 1 nearest neighbor Heisenberg antiferromagnet is discussed 
as it is tuned from a two-dimensional square lattice, of lattice spacing a, towards a spin chain by 
varying the width L y of a semi-infinite strip L x x L y . For integer spins and arbitrary L y , and for 
half integer spins with L y /a an arbitrary even integer, explicit analytical expressions for the zero 
temperature correlation length and the spin gap are given. For half integer spins and L y /a an odd 
inetger, it is shown that the c = 1 behavior of the SU (2)i WZW fixed point is squeezed out as the 
width L y — > oo; here c is the conformal charge. The results specialized to S = 1/2 are relevant to 
spin-ladder systems. 



PACS: 75.10.Jm, 75.10.-b 

One-dimensional quantum antiferromagnets have 
many unusual properties. For example, nearest-neighbor 
Heisenberg spin chains with half integer spins are gap- 
less, but those with integer spins are generically gapful 
11. Many properties, such as these, can be understood 
from the correspondence of the spin chain to a (1 + 1)- 
dimensional 0(3) non-linear cr-model, augmented by a 
term in the action that contributes a phase factor e l8 ® 
to the path integral, where = 2nS ; S is the spin, 
and Q is an integer winding number. For integer spins, 
6 = mod (2ir), the phase factor is unity, but for half- 
integer spins, = 7T mod (2ir), it is (— 1)^. This leads to 
the crucial difference between the integer and the half- 
integer spins. Although the existence of the gap in the 
model with = can be seen in a variety of ways [Q, 
the model with 9 = n is more subtle ||. For 9 = n, the 
short distance behavior is dominated by two weakly in- 
teracting goldstone modes. These are correctly described 
by the perturbative renormalization group analysis that 
is impervious to the existence of the #-term. In the lan- 
guage of conformal field theory, the system is in the prox- 
imity of an infrared unstable fixed point corresponding 
to conformal charge c = 2. At short distances, there is 
no distinction between integer and half integer spins, and 
the system appears gapful. However, on longer scales, the 
9 = 7r model flows to k = 1 SU(2) Wess-Zumino-Witten 
model (SU(2)i WZW), corresponding to a c = 1 mass- 
less theory. Indeed, all critical theories in two dimensions 
must belong to a conformally invariant fixed point. 

The two-dimensional square-lattice, nearest neighbor 
Heisenberg antiferromagnet is entirely different. It is rig- 
orously known that the ground state is Neel ordered for 
S > 1 Q, but no such proofs exist for S — 1/2. Nonethe- 
less, the numerical evidence for an ordered ground state 
for S = 1/2 is strong ||. Moreover, the assumption of an 
ordered ground state has yielded predictions B that are 
confirmed in neutron scattering experiments pj. Hence- 
forth, I shall assume that this is also a solved problem, 
and the correct low energy theory is a (2 + l)-dimensional 



0(3) non-linear cr-model, which is in the broken symme- 
try phase in its ground state || . (The possible existence 
of topological terms were considered and discarded by 
a number of authors ||].) The corresponding elementary 
excitations are weakly interacting Goldstone modes, that 
is spin waves. This low-energy, long wavelength model is 
essentially geometrical and is almost entirely determined 
by symmetry || regardless of whether or not the mag- 
nitude of the spin is large. The two needed phenomeno- 
logical constants are the spin wave velocity and the spin 
stiffness constant. In principle, experiments can deter- 
mine these constants and there is no need to rely on a 
presumed large-S* expansion. 

It is an intriguing question to ask how a two dimen- 
sional system would evolve if we began with a strip, 
L x x L y , and continuously tuned the system by vary- 
ing the width L yi with L x kept infinitely large. Would 
it approach the one-dimensional limit, and, if so, how 
would we recover the sensitivity to the topological an- 
gle 91 This would be of only theoretical interest, albeit 
considerable interest, if it were not for experiments on 
spin ladders (l(J in which S — 1/2 systems of varying 
width are explored. The purpose of the present paper is 
an attempt to clarify this crossover and to determine the 
evolution of the excitation spectrum. 

It has been argued [ pd|Jl2| that spin ladders correspond 
to an effective (1 + l)-dimensional 0(3) non-linear cr- 
model with the ^-parameter given by 9 = 2imiS, where 
ri; is the number of legs. Thus, for S = 1/2, the system is 
gapful for even-leg ladders and gapless for odd-leg ladders 
in accord with experiments [fl0|| . There are complimen- 
tary theoretical and numerical approaches to the ladder 
systems that are outside the scope of the reasoning in 
this paper [^3|. However, the crossover problem stated 
above has not been fully elucidated, although it was an- 
ticipated pl| that, when approached along the sequence 
of even-leg ladders, the gap must collapse exponentially 
with the increasing width of the system. In the present 
paper, I shall show precisely how this happens and de- 
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rive a formula that can be checked. At first sight, the 
approach to the two-dimensional limit along the odd se- 
quence appears to be simple, because they are gapless. 
However, this is not so because the two-dimensional prob- 
lem, which is insensitive to topology, is described by two 
Goldstone modes; it cannot be a straightforward exten- 
sion of the c = 1 fixed point of the SU{2)\ WZW model. 

The Euclidean action of the O(N) quantum non-linear 
cr-model, in which one of the dimensions is singled out 
and is of finite extent L, is 



1 = 2%' dT 



jd-l 



x i dx\ 





(1) 



where the index p runs over all the spatial dimensions, 
1 through d. The extent of the imaginary time dimen- 
sion, (3%, tends to infinity as the temperature T — 1/ksP 
tends to zero. We shall impose a periodic boundary con- 
dition along the direction 1; the remaining spatial di- 
rections will be assumed to be infinite in extent. The 
staggered order parameter field of the antiferromagnet, 
f2, is an A-component unit vector field, which is a func- 
tion of (r, x\, X2, ■ ■ ■ , Xd)] the spin-S 1 antiferromagnet cor- 
responds to A = 3. The parameter p° s is the bare spin 
stiffness constant at the spatial cutoff, A -1 , of the model, 
and the parameter c is the spin wave velocity on the same 
scale. I shall focus on the zero temperature behavior and 
report results of a "Lorentz" invariant analysis; therefore, 
the spin wave velocity will not renormalizc. 

The action at T = is interesting. The extents in 
all directions, except x\, are infinite; along ii we have 
a periodic boundary condition. The physical problem 
at T = is, therefore, isomorphic to a problem at fi- 
nite "temperature" , where the temperature-like variable 
is el = Tic/ L. With proper identifications of the parame- 
ters, it is identical to that solved in Ref. Let us define 
two dimcnsionlcss bare coupling constants: 
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The energy like parameter e l plays the role of the di- 
mensionless temperature-like coupling constant in Ref. 

H , and go is the same as that defined previously. 

The renormalization group equations can be simply 
read off from Ref. ||. The crossover phase diagram, in 
d = 2, constructed from these equations is shown in Fig. 

I, merely for orientation. The three distinct regions had 
to be renamed, as the present analysis corresponds to 
T = 0. The regions previously named "renormalized clas- 
sical" , "quantum critical" , and "quantum disordered" 
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FIG. 1. The crossover phase diagram in d — 2. 

are now renamed to be "dimensionally reduced" , "criti- 
cal spin liquid", and "gapped spin liquid", respectively. 
With simple transcriptions, the physical pictures of the 
crossover boundaries are the same as before. The anal- 
yses of the gapped and the critical spin liquid regimes 
are exceedingly complex and are beyond the scope of the 
present paper. 

The region of the phase diagram for which we can make 
precise predictions is the dimensionally reduced region. 
In this region, the system is in the Neel ordered state 
when L — oo, or el — 0. When L ^ oo, the system is 
equivalent to a dimensionally reduced effective (1 + 1)- 
dimensional model with no long range order. This can 
be seen from the renormalization group equations . At 
first, with increasing length scale, g rapidly decreases, 
and El increases slowly, that is, the system appears more 
and more ordered. Subsequently, the growth of g slows 
down, but El increases more rapidly, thereby breaking 
up the order at longer length scales, resulting in the re- 
duction from (2 + 1) to (1 + 1) dimensions. The effective 
coupling constant to be used in the (1 + l)-dimensional 
model is easily calculated to be 
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where p s is the fully renormalized macroscopic spin stiff- 
ness constant at T = of the square-lattice spin-S 1 an- 
tiferromagnetic Heisenberg model (L x = oo, L y = oo) 
|]l4j| . This definition can be made more explicit, if we 
recall that p s = JS 2 Z Ps , and Tic — 2^/2JSaZ c , where 
J is the exchange constant, a is the lattice spacing; Z Ps 
and Z c are the renormalization factors |Q . We can now 
write 



£cB -s 



V2Z C J a ttS 



(4) 



Therefore, for large (L/a), the input bare coupling con- 
stant to the efffective (l + l)-dimensional model is greatly 
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reduced from its value (2/5) of a spin chain. This is due 
to increased order at short distances, concommitant of 
the quasi two-dimensional nature of the model with fi- 
nite width L. 

For integer spins, and for even-legged ladders with half 
integer spins, the description of the dimensional crossover 
is conceptually complete. The system is massive and its 
mass gap should be calculated with the effective cou- 
pling derived above. I shall make more precise predic- 
tions later. The case of odd-legged ladders with half 
integer spins requires further clarification. Because all 
possible topological terms were dropped in the (2 + 1)- 
dimensional model, the masslessness of the dimensionally 
reduced system could not be recovered. Even in the pres- 
ence of sufficiently strong local Neel order, the proof of 
the nonexistence of the topological term in d = 2 is cor- 
rect strictly when the number of spins along both L x and 
L y are even. For odd number of chains along L y , but 
L x = oo, a topological term 2ttiSQ remains ||. If wc 
wish, we can rewrite it as 2i:iSQni by realizing that the 
topological angle is only defined modulo 27r. It has been 
never fully explained why this odd-chain case is physi- 
cally irrelevant; in fact, it is not, as we shall see. Imagine 
that we include such a term in our action fil| . For finite 
L, this should, in principle, render the model massless in 
the sense of c ~ 1, not in the sense of a Goldstone phase. 
However, the c = 1 behavior will be difficult to see when 
L/a is large. The reason is that the effective coupling 
constant derived in Eq. (4) will be very small, and the 
perturbative renormalization group, which is impervious 
to the topological term, will be valid up to very long 
distances, until the coupling constant becomes of order 
unity for the system to crossover to c = 1 SU (2)i WZW 
model. Recall that the 9 = tt model in (1 + 1) dimensions 
is massless but not conformally invariant; it has a non- 
trivial /3-function and an associated mass gap. Thus, the 
region in which the c = 1 behavior is seen is squeezed 
out as L — > oo, and all we see is the Goldstone phase; 
we shall see that the mass gap vanishes exponentially as 
L — > oo. Conversely, when (L/a) is of order unity, the 
c = 1 feature should be visible. 

Using Ref. II , it is possible to write down by inspection 
the expression for the correlation length, £, in our 0(3) 
model defined on a strip for which L x = oo, L y = L. 
It is important to note that the dimensionally reduced 
effective (1 + l)-dimensional model has "Lorentz" invari- 
ance; integrating out the Ly-modes does not destroy the 
proportionality between the imaginary time and the L x 
directions. There is, therefore, one and only one correla- 
tion length . The result for £ is 
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From strong coupling simulations, the quantity 27rC was 
estimated in Ref. [|| to be between 0.01 and 0.013. This 
makes the overall numerical prefactor to be between 0.27 
and 0.35. Since then an asymptotically exact expression 
has been derived and the exact prefactor is known 
to be e/8 » 0.34. In addition, the constant which was 
previously known to be only of order unity is determined 
to be A = 1/2. 

It is interesting to rewrite Eq. (5) in terms of the 
Josephson correlation length of the L = oo (2 + 1)- 
dimensional 0(3) non-linear cr-model p|, which is given 
by £j = (he/ ps). This length separates the short dis- 
tance critical behavior from the long distance Goldstone 
behavior. In terms of £j, the correlation length takes the 
simple finite-size scaling form: 

e = Kfe) e ^ L/t ' i 1 s {c,/27tl) + ° { ° /2nL) ' 

(6) 

For spin-5 square- lattice Heisenberg antiferromagnet, £j 
is given by [[l5| 



( 2^2Z C 



(7) 



(5) 



Because of "Lorentz" invariance, the spin gap, A, is 
simply A = hc/t;. Note again that the spin wave velocity 
does not renormalizc and this relation does not have any 
corrections (cf. below). Specializing to 5 = 1/2, we get 

(£/a)i = 0.499e a682(i/Q) [l- 0.734(a/L)], (8) 
(A/J), = 3.347e- a682(i/a) [1 - 0.734(a/L)]" 1 . (9) 

The field theory analysis presumes that a continuum the- 
ory is applicable in both directions, and therefore the 
expressions in Eqs. (8) and (9) cannot be accurate for 
(L/a) ~ 1. Moreover, these expressions, obtained with 
a periodic boundary condition, are likely to be differ- 
ent from those obtained from other boundary conditions, 
such as the open boundary condition, especially when 
(L/a) ~ 1. Nonetheless, if wc take L/a — 4 correspond- 
ing to four-leg ladders, we get (A/J)i = 0.268, and 
(£/a)i = 6.23. For L/a = 6, we get (A/J)i = 0.064, 
and (£/a)i = 26.2. 

To compare, numerical results, with open boundary 
condition along L y , are available for L/a = 4 and 6. 
According to Ref. [[ITJ, we have, for L/a = 4, (A/J)i = 
0.190, and (£/a)i =5-6. According to Ref. 0,'we 
have, for L/a = 4, (A/J)i = 0.160, and (£/a)i = 10.3. 
From the same work, we have, for L/a = 6, (A/J)i = 
0.055. The value for the correlation length for L/a = 6 
is not given explicitly. However, a simple extrapolation 
yields a value close to 30. From Ref. pl[ , we have, for 
L/a = 4, (A/J) i = 0.17; for L/a = 6, (A/J) i = 0.05. 
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Note that the gaps should be inversely proportional to 
the correlation lengths, where the proportionality con- 
stant does not renormalize due to "Lorentz" invariance. 
We require "Lorentz" invariance only in the L x -t plane. 
This has little to do with the fact that the width along 
L y is finite and equal to L. This proportionality is auto- 
matically satisfied for the analytical expressions given in 
this paper and should be a good check on the numerical 
work. 

Based on the simple observation that at T — the 
spin-5 square-lattice Heisenberg model of finite width 
can be mapped on to a (2 + l)-dimensional 0(3) non- 
linear cr-model with a finite dimension, I have provided a 
theory for crossover in spin ladders. Explicit analytical 
expressions for the correlation length and the spin gap 
were obtained by transcribing the results in Ref. || . The 
agreement with numerical calculations is very good con- 
sidering the difference in the boundary conditions em- 
ployed. The analytical expressions show precisely that 
the crossover to the two-dimensional limit is approached 
exponentially as the the number of legs in the ladder sys- 
tem is increased. The extension of the present theory to 
anisotropic coupling and to finite temperature should be 
straightforward. I hope to return to these extensions in 
the future. 

This work was conceived and carried out at the Aspen 
Center for Physics. It was supported by a grant from the 
National Science Foundation, Grant No. DMR-9531575. 
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